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ABSTRACT 

We have examined the electroweak radiative corrections in the LEP precision 
data in view of the new measurements of Mw and rut as well as the recent 
progress in the higher order radiative corrections. From the minimal x^-^i to 
the experimental Z-decay parameters (with the aid of a modified ZFITTER 
program), we predict that My/ = 80.29(4) (2) GeV where the first error is due to 
the uncertainty in the fitted mt for a fixed mn and the second error comes from 
the tuh in the range of 60 — 1000 GeV, which is to be compared with the current 
world average Mw = 80.23(18) GeV. The current world average value of Myy 
and the 1994 LEP data definitely favor nonvanishing electroweak radiative 
corrections and are consistent with a heavy mt as measured by the recent CDF 
report but with a heavy Higgs scalar of about 400 GeV within the context of 
the minimal standard model. The sensitivity of and the errors in the best fit 
solutions due to the uncertainties in the gluonic coupling as{Mz) and a{Mz) 
are also studied carefully. In addition we discuss how the future precision 
measurements of Mw can provide a decisive test for the standard model with 
radiative corrections and give a profound implication for the measurement of 
t-quark and Higgs masses. 



^Supported in part by the USDOE contract DE-FG02-91ER40688-Task A. 

^ Supported in part by the Basic Science Research Institute program, Ministry of Education, 1994, 
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Much interests have been paid in recent years to the electroweak radiative corrections 
(EWRC) and precision tests of the standard model thanks to the accurate data obtained 
at LEP [1-7]. There have been numerous articles published on the subject as has been 
documented in [7-9]. The LEP data are generally regarded as the support for the standard 
model and as the evidence of the nonvanishing EWRC [9]. 

However there have been several important developments since last year which war- 
rant a new motivation to repeat the precision tests of the standard model. Some of the 
experimental advances are (1) the new CDF value for [10], (2) the improved LEP 
data [11], and (3) the CDF and DO reports 0on nif [12], while there have been also some 
progress on the higher order corrections, in particular the dominant two-loop terms of or- 
der a'^mj [14], the QCD-electroweak mixed diagrams [14,15] and higher order corrections 
to the QCD factor in the Z- decay width [16]. We would like to present the results of the 
new fit to the updated 1994 data with the aid of the appropriately modified ZFITTER 
program to incorporate these new theoretical developments. We examine the errors in 
the best fit solutions due to the uncertainties in the strong coupling constant as{Mz) 
and also in a{Mz)- In the analysis we determine Mw self-consistently from the W-mass 
relation that includes EWRC for the value of mt covering experimental range and fit the 
LEP data, and show how stable the predicted Mw is regardless the exact value of mn 
in the interesting range of 60 — 1000 GeV. Though the sensitivity of the EWRC to the 
exact value of Mw in the standard model has been studied based on the ly-mass formula 
[17], the effect of the self-consistency satisfied by M\y through the mass relation with 
EWRC as well as that of the uncertainties in as{Mz) and a{Mz) to the precision tests 
of the electroweak data has not been fully examined and understood. For this reason, 
we would like to critically examine in this paper the sensitivity of the precision tests and 
the rut — niH correlation to the requirement of consistency in the needed EWRC for a 
set of mass values as well as the errors in as(M^) and a{Mz)- In particular, the results 
of the minimal x^"fit show that the CDF value = 174 GeV can be consistent with 
a best fit solution but with a Higgs mass about mn = 400 GeV if as{Mz) = 0.123 and 
a-^{Mz) = 128.87 and predict Mw = 80.29(4) (2) GeV where the first error is due to the 
uncertainty in the fitted for a fixed mu and the second error comes from the tuh in the 
range 60 — 1000 GeV. However mt and mn can easily be shifted, due to the uncertainty 
in the gluonic coupling Aa^ = ±0.006, by as much as 5 GeV and 125 GeV respectively, 
while the corresponding shift in Mw turns out to be about 30 MeV. The possible shifts 
in mt, mu and Mw due to the error in the gauge coupling Aa;~^(M^) = ±0.12 can be as 
much as 6 GeV, 160 GeV and 20 MeV respectively. 

3The most recent values [13] are mt = 176±8±10 GeV (CDF) and mt = 199l:2?±22 GeV (DO). 
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In addition, we reexamined the claim made by Novikov, Okun, and Vysotsky [18] 
that the 1993 data from LEP on the electroweak parameters as defined in the standard 
model could be explained by the QED Born approximation (QBA) in which a{Mz) is used 
instead of a(0) in the tree approximation along with the corresponding redefinition of the 
weak mixing angle sin^ 9 instead of sin^ 9w- In particular the so-called QBA predictions 
were claimed to be within la accuracy of all electroweak precision measurements made 
at LEP in 1993. In order to examine the intriguing claim made in Ref.[18], we considered 
the case of the QBA by consistently neglecting the terms of non-photonic origin in the full 
EWRC. The full EWRC is calculated with the aid of a modified ZFITTER program [19] 
that uses an improved QCD correction factor [16], includes the dominant two-loop and 
QCD-electroweak terms [14] and makes the minimal fit to the data. We have found 
[21] that the situation is very sensitive to the value of Mw and that the QBA fit to the 
1993 data gives statistically comparable or better and therefore the existence of the 
genuine quantum effect due to full EWRC might not have been evident in the 1993 data. 
We note that the old M]^ has somewhat lower central value and larger error. The new 
results of the minimal fit to the updated 1994 data and the predicted M\y compared 
to the new world average Mw show a clear effect due to the EWRC. 

It is well known that the charge renormalization in the conventional QED fixes 
the counter term by the renormalized vacuum polarization n'''(0) and one can evaluate 
fPlq"^) = TP^{q^)/q^ from the photon self energy il'^^{q^), for example, by the dimensional 
regularization method. This gives at q^ = M| the total fermionic contribution oi uif < 
Mz to the real part RetViyM"^) = —0.0596(9), which includes both the lepton and quark 
parts [22]. Here, the quark contribution to Rell''{q'^) is the hadronic one which can 
be directly evaluated by dispersion integral over the measured cross section of e~^e^ 
hadrons. Then, we get a{Mz) = 1/128.87(12) in the on-shell scheme if the hyperfine 
structure constant a = e^/47r = 1/137.0359895(61) is used, which we will use in this 
paper. The error in a{Mz) is due to the uncertainty in hadronic contribution. This is 
obviously a source of the uncertainty in the best fit solutions. 

The electroweak parameters are evaluated numerically with the hyperfine structure 

constant a, the four-fermion coupling constant of /i-decay, = 1.16639(2) x lO^^GeV"^, 

and Z-mass Mz = 91.1888(44) in the 1994 data fit. Numerical estimate of the full EWRC 

requires the mass values of the leptons, quarks, and Higgs scalar besides these quantities. 

While Z-mass is known to an incredible accuracy from the LEP experiments largely 

due to the resonant depolarization method, the situation with respect to the l^-mass is 

''We note that the higher QCD effects (aa^mj order) are yet to be unanimously agreed by the experts 
as evidenced by the discussions in [15,20]. 
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Mw (GeV) 


Ar 


sin^ 9w 


1 


79.91 


0.0623 


0.2321 


2 


80.22 


0.0448 


0.2261 


3 


80.23 


0.0443 


0.2263 


1 


80.38 


0.0355 


0.2231 



Table 1: Dependence of the radiative correction Ar on the values of Mw ■ 



desired to be improved, i.e., Mw = 80.22(26) GeV [23] and 80.23(18) GeV[10] from the 
world average values vs. the old and new CDF measurements Mw — 79.91(39) GeV 
[24] and Mw = 80.38(23) GeV[10]. The minimal x^-fit to the LEP data will at best 
give the mt — mn correlations as shown in Fig. 1(a) and Fig. 1(b). We will show how 
the best fit solutions are chosen out of the solution set for (mt,mij) and determine Mw 
self-consistently from the 1^-mass relation that includes EWRC and how they depend on 
the uncertainties in as{Mz) and a{Mz) so that how precisely the standard model can 
be tested at the moment, let alone the predictions of the rrit — mn correlation from the 
experimental M^r. 

One has, in the standard model, the on-shell relation sin^ 9w = 1 — while the 
four-fermion coupling constant can be written as 



z 



na ( ^ 




so that Ar, representing the radiative corrections, is given by 

/37.2802\2 1 

We note from Tabic 1 that the radiative correction Ar is very sensitive to the value of 
Mw- Merc change in Mw by 0.59% results as much as a 43% change in Ar. Theoretically, 
the radiative correction parameter Ar within the standard model can be written as [25] 

1 - Ar = (1 - Aa) • (1 + cot^ %Ap) - Ar^em, (3) 

where Ap contains one loop and the leading 2-loop irreducible weak and QCD corrections. 
Main contribution to Ap = 1 — p^^ is from the heavy t-quark through the mass rcnormal- 
izations of weak gauge bosons W and Z, while there is a part in (Ar)rem containing also the 
t-quark and Higgs scalar contributions. Note that the so-called QBA to Ar is defined by 
keeping only the photon vacuum polarization contribution, Aa = — i?efl'^(M|) = 0.0596. 
We see from Table 1 that Aa is numerically the dominant component of the radiative 
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corrections. In particular for the old CDF Mw, Aa is already within 4.4% of the needed 
Ar and is close enough to be within the experimental uncertainty. However, with the 
current world average value Mw — 80.23 GeV, Aa differs by 32% from the needed Ar 
that has to be accounted for by the weak interaction corrections [26] . 

Note that precise determination of the on-shell value of sin^ 9w can also constrain 
the needed radiative correction and the value of Mw 

We have searched for the minimal x^-fits to both 1993 and earher 1994-data of the Z- 
decay parameters measured at LEP by using the modified ZFITTER program [21] . Within 
the framework of the standard model in which G^, a and Mz are taken as input, one can 
predict Mw from (2) in a self-consistent manner. Starting with the given masses of the 
quarks, leptons, gauge bosons and Higgs scalar, as well as given gluonic couphng as{Mz) 
and gauge coupling a{Mz), Ar is calculated from (3) by including up to the dominant 
two-loop and QCD-electroweak mixed terms and then is used to determine Mw from the 
right hand side of (2). With this new Mw, Ar is calculated again to determine another 
new Mw- This iteration process is repeated until Ar converges to within O(10~^). The 
final output Mw from this iteration procedure is the self-consistent solution of (2) for Mw 
with the starting set of rrit, uih etc. Upon varying rrit, this procedure will give the rrit—Mw 
correlation for all other parameters including niH fixed. The family of the curves in Fig. 
2 represents such correlations for different fixed values of niH- The rrit ~ "^i? correlation 
shown in Fig. 1(a) and Fig. 1(b) is not restrictive enough to discriminate an interesting 
region of rrit and run when compared to the experimental range of rrit- Also the Mw — frit 
correlation shown in Fig. 2 does not discriminate mu as long as it is heavier than 100 
GeV given the current experimental situation of {rrit, ^w)- Fig- 3 shows the result of the 
self-consistent procedure for Mw — mn correlation when rrit is restricted to the CDF value 
174 ± 16 GeV on which the predicted Mw range is also indicated. The error band ±0.04 
GeV in the predicted Mw is due to the uncertainty in the fitted rrit for a fixed mn as one 
can deduce from Table 2. We then calculate the eleven Z-decay parameters, as chosen in 
Table 2, for the parameter sets {rrit, t^h) that determine Mw from (2) and search for the 
minimal x^-fit solution to the experimental Z-decay parameters. This procedure selects 
the Best. fit curve in Fig. 1 and o points on each curve in Fig. 2. Including the error due 
to varying rriH in the range, we predict Mw — 80.29(4) (2) GeV, which is shghtly larger 
than the central value of the current world average 80.23 GeV. Note that the predicted 
Mw constrains to he 230 GeV and 830 GeV for rrit — 174 ± 16 GeV as one can see 
from Fig. 3. 

The results of the best x^-^ts to the updated 1994 data are given for the four 
sets of {mt,mH) in Table 2. The Z-decay parameters are calculated with the gluonic 
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coupling constant in the range as{Mz) — 0.123(6) in the QCD correction factor Rqqj) — 

1 + 1.05^ + 0.9(±0.1) {^^y - 13.0 for light quarks(u,d,c,s) [27] and in Rqcj) = 

1 + ci(mfe)^ + C2{mt,mt) (^^^ + cs{mt,mb) , the mb and rrit mass dependent one 
for b quarks [16]. The partial width for Z — > // is given by 

T/ = ^||/3flQEDC/flQCD(M|) + [aff] x + 2^) - (4) 

where (3 = ^1 — Arnj/M^, -^QED = ^ + If '^f ^^id the color factor cj = 3 for quarks 
and 1 for Icptons. Here the renormalizcd vector and axial-vector couplings arc defined 
by aj = ypj2aj = ^J~f^2ll and = aj[l — 4|Qj| s\v? O^k^] in terms of the familiar 
notations [19,28]. Note that Aa is included in the couplings through siv? 6^ via (1) and 
(3) and all other non-photonic loop corrections are grouped in and as in [19,29] 
including the dominant two-loop and QCD-electroweak terms. Note that the QED loop 
corrections can unambiguously be separated from the electroweak loops in the case of 
neutral current interactions [28]. Thus the case of the QBA can be achieved simply by 
setting and to 1 in the vector and axial-vector couplings. 

Numerical results for the best x^-fits to the 1993 LEP experimental parameters of Z- 
decay for Mw = 79.91(39) GeV [24] and Mw = 80.22(26) GeV [23] showed [21] generally 
small contributions of the weak corrections and in particular that the QBA was close to 
the experimental values within the uncertainty of the measurements, i.e., within 2a. The 
near absence of the weak interaction contributions to the radiative corrections for the 1993 
data is more impressive for Mw = 79.91(39) GeV than for Mw = 80.22(26) GeV. This 
is mainly because QBA gives Mw — 79.95 GeV compared to Mw = 80.10(1) GeV from 
the full EWRC for mu in the range of 60 — 1000 GeV and a larger error in the former as 
observed in [18,21]. At closer examination, however, the QBA in this case over-estimates 
the radiative corrections and the full EWRC fair better; for Mw = 80.22(26) GeV and 
mjj = 60 - 1000 GeV, Ar = 0.0596 in QBA and 0.0498 - 0.0505 in the full EWRC to 
be compared to the required value 0.0448. Also the global fit to the 1993 data with two 
variables rrit and ttih in the range 60 — 1000 GeV show [21] that the best fits can be 
achieved by nit = 139(17) GeV wth a stable output Mw = 80.12(1) GeV when the full 
EWRC are taken into account. 

The situation with the minimal x^-fits to the updated 1994 LEP data and with 
theoretically determined Mw is significantly different from the case of the 1993 data as 
one can see from Table 2. Not only there is clear evidence of the full EWRC in each of 
the eleven Z-parameters but also the best fit solutions to the 1994 data show a stable 
output Mw = 80.29(2) GeV for mn in the range of 60 - 1000 GeV. In particular the 
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Table 2: Numerical results including full EWRC for eleven experimental parameters of 
the Z-decay and Mw- Each pair of mt and rriH represents the case of the best x^' fit to 
the 1994 LEP data for a^(Mz) = 0.123(6) and a'^Mz) = 128.87(12). The numbers in 
represent the errors due to Aas{Mz) = ±0.006 and Aa~^{Mz) — ±0.12 respectively, 
sin^^gj^^"" is from the measurement of < Qfb >■ For the case of Born approxiamtion, 
the errors are due to Aa~^{Mz) only. 
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QBA gives distinctively inferior x^(=19.5/ll) for the 1994 data. Also the CDF rrit = 174 
GeV is a possible output solution with a mn about 400 GeV among the many possible 
combinations of {mt,mH) given by the 'Best. fit' curve in Fig. 1(a) and Fig. 1(b). As 
shown in Table 2, the Best.fit solutions can have errors due to the uncertainty in as{Mz) : 
nit and Mw may be shifted by as much as ±5 GeV and ±30 MeV respectively because of 
Aag — ±0.006. The error range of the Best.fit solutions is indicated by the curves A and 
B in Fig. 1(a). There are additional comparable errors due to the uncertainty in a{Mz) 
as shown in Table 2 and Fig. 1(b) : Aa~^{Mz) — ±0.12 can cause another ±6 GeV and 
±20 MeV respectively in rrit and Mw. In general the x^-value tends to prefer the lower 
rrit and accordingly smaller niH, though any pair of {rrit, ttih) on the Best.fit curve in Fig. 
1(a) and Fig. 1(b) is statistically comparable to each other. In particular the best global 
fits to the updated 1994 data give m* = 155 - 187 GeV for mn = 100 - 1000 GeV. Most 
of the Z-parameters are stable irrespectively to the uncertainties due to Aag and Aa and 
in excellent agreement with the data except R{riji/rfiad)- Even with the mass dependent 
QCD factor, there is still about 2.4 a deviation in i?(rj5/r^ad) from the experiments 
irrespectively to the uncertainties in as{Mz) [15]. Most of the contributions are from 
R{r^l/rhad) and to a lesser degree from Rircc/I'had) and A^g. Fig. 2 shows how Mw 
changes with rrit for fixed mn from the consistency of the full EWRC, on which the new 
world average Mw and CDF rrit along with the Best.fit solutions (depicted by o points) 
are also shown. The central values of the world average Mw and CDF rrit are consistent 
with a Higgs scalar mass about 1000 GeV, though mn — 100 GeV is within 1 a because 
of large errors in the data. Clearly a better precision measurement of Mw is desired to 
distinguish different mn- For example, a change of mn by 200 GeV, i.e., from 400 GeV 
to 200 GeV at rrit = 174 GeV, results a change of 50 MeV in Mw, i.e., from 80.30 GeV to 
80.35 GeV, as one can see from Fig. 2 and Fig. 3 . This in turn will require a precision 
of 11 GeV or better in rrit from the Best.fit curve in Fig. 1(a), which is consistent with 
the most statistical error improvement that may be achieved at the Fermilab Tevatron. 
Present precisions in the data entail a theoretical uncertainty of about 36 MeV in Mw 
which is about the overall error improvement expected at LEP-200. 

We have examined the results of the minimal x^-fits to the precision measurements of 
the Z-decay parameters at LEP with the aid of a modified ZFITTER program containing 
the full one-loop and dominant two-loop EWRC. While the result of QBA might appear 
to be in agreement with the 1993 data within 2a level of accuracy [21,30], the new world 
average value of Mw and updated 1994 LEP data definitely disfavor the QBA and support 
for the evidence of the nonvanishing weak correction. Furthermore, the CDF rrit is a 
solution of the minimal x^-fi^s to the 1994 data with a Higgs scalar mass about 400 GeV. 
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However this rrit value can be shifted by as much as 8 GeV due to the overall uncertainties 
in as{Mz) and a{Mz) for the moment and accordingly mn ranging 250—690 GeV. Further 
precision measurement of can provide a real test of the standard model as it will give 
a tight constraint for the needed amount of the EWRC and can provide a profound 
implication for the mass of t-quark and Higgs scalar. If Mw is determined to within a 40 
MeV uncertainty, Ar within the context of the standard model will be tightly constrained 
to distinguish the radiative corrections and the x^-fit to the Z-decay data with the 1994 
accuracy can discriminate the mass range of the t-quark and Higgs scalar within 8 GeV 
and 200 GeV respectively, providing a crucial test for and even the need of new physics 
beyond the standard model. Finally if Mw is determined to be larger than 80.31 GeV with 
better than a 20 MeV accuracy by the future precision measurements (perhaps reachable 
at LHC), this would be a definite sign for new physics beyond the standard model. 

Note added: Most recent CDF measurement is Mw = 80.41(18) GeV [31]. We note 
that the theoretical prediction of Mw, 80.29(4) (2) GeV, is about 0.67(7 below the central 
value of the new experimental Mw 
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Figure Captions 



Fig. 1(a) : The mass ranges of rut and mn from the minimal x^-fit to the updated 1994 
LEP data for as{Mz) = 0.123(6) and a-^{Mz) = 128.87(12). The range of the 
Best. fit solutions (see the text) due to the error Actg = ±0.006 is also indicated by 
the curves A and B. 

Fig. 1(b) : The mass ranges of rrit and itlh from the minimal x^"fit to the updated 
1994 LEP data for ^^(M^) = 0.123(6) and a-\Mz) = 128.87(12). The range of 
the Best. fit solutions due to the error Aa'^{Mz) = ±0.12 is also indicated by the 
curves a and b. 

Fig. 2 : Mw versus rrit for fixed values of niH from the full radiative correction in the 
standard model. The case of the minimal x^-fit to the updated 1994 LEP data are 
indicated by O. 

Fig. 3 : M]y versus rriH for rrit = 174± 16 GeV : The middle region bounded by the solid 
lines represent the predicted theoretical value of My[/. Notice that as rrit becomes 
heavier than 174 GeV the lower bound of thh increases above 300 GeV from the 
predicted Mw- 
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